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l_i^ , Employing different statistical ensembles may lead to qualitatively different results concerning 

averages of physical observables on the mesoscopic scale. Here we discuss differences between the 
canonical and the grandcanonical ensembles due to both quenched disorder and thermodynamical 
effects. We show how these differences are related to spectral correlations of the system at hand, and 

1.^ ' evaluate the conditions (temperature, system's size) when the thermodynamic limit is achieved. We 

demonstrate our approach by evaluating the heat capacity, persistent currents and the occupation 
probability of single electron states, employing a systematic diagrammatic approach. 
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I. INTRODUCTION 

Standard textbooks on equilibrium statistical mechanics often assert that the various statistical ensembles in stan- 

q_) ' dard use (e.g., the canonical ensemble (CE) and the grandcanonical ensemble (GCE) are equivalent. Landau and 

Ch Lifshits [1] assert hat ". . . all three distributions, the microcanonical and the two forms of the Gibbs distribution, are 

C^ in principle suitable for determining the thermodynamic properties of the body. The only differnce from this point 

Ch ■ of view lies in the degree of mathematical convenience" . Huang [2] claims that ". . . the grand canonical ensemble is 

^ trivially equivalent to the canonical ensemble for N particles" , and there are many more examples. These assertions 

O should be understood as referring to average thermodynamic quantities (as distinct from fluctuations; it is evident 

'— '^ that by the very nature of the respective ensembles, fluctuations in the number of particles, for example, are very 

^~~l ■ different between the CE nd the GCE). Moreover, it is assumed that the body under consideration is macroscopic 

OO ' (i-G., the "thermodynamical limit" has been taken). In this limit the relative fluctuations tend to zero. 

p^ ' It is important to stress, though, that thermodynamic (or statistical-mechanical) quantities may be perfectly defined 

i^ [ even for a system with few degrees of freedom - provided they are in equilibrium with an infinite reservoir (heat 

r ■ I bath, particle reservoir, etc.). A well-known example is the Fermi-Dirac distribution function which is obtained by 

considering a single-level system coupled (weakly) to an electron bath at temperature T and chemical potential /i. In 

the limit of small systems fluctuations are large (relative to mean values). As we have stressed above, their nature 

_JL ' is very different between the different ensembles. More interestingly, average thermodynamic quantities too may 

G I be different between the CE and the GCE. Most strikingly, qualitative differences between the two ensembles may 

O ■ arise when quantum coherence comes into play, on mesoscopic scales and downwards. Novel effects may appear. 
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which are quite remote from our common wisdom and everyday intuition as applied to macroscopic systems in the 
thermodynamic limit. The thermodynamic limit is attained as the system's size tends to infinity and/or when the 
temperature is made sufficiently high. 

It is evidently interesting to identify the scale of temperature (in conjunction with the system's size) at which 
the non-trivial differences (i.e., differences in mean values, possibly qualitative) between the CE and the GCE are 
suppressed. Furthermore, we shall note that on finite length or temperature scales there are certain differences 
between thermodynamic relations and statistical-mechanical ones (therefore these two different frameworks may lead 
to different predictions, although, as was noted above, both statistical mechanics and thermodynamics are well- 

1 



defined even in ttie context of finite systems. Consequently there are cliaracteristic temperature or length scales which 
mark the suppression of the differences between statistical mechanics and thermodynamics, and the crossover to the 
"thermodynamic limit" . 

It should be stressed that in most practical instances the fabrication of samples is not a perfectly reproducible 
process. Commonly there are sample-to-sample variations in sample size and shape, in the mcsoscopic impurity 
configuration and the location defects (we shall refer to all these factors as "disorder" ) . Evidently it is possible to ask 
questions pertaining to sample specific obscrvablcs, and to evaluate the magnitude of sample-to-sample fluctuations 
in these obscrvablcs. But it is also of interest to consider ensemble averaged quantities, where averaging over disorder 
is included (the concept of ensemble will be elucidated below) . Such ensemble averaged quantities will be at the focus 
of the present study. 

The purpose of this paper is twofold. We shall first review and expand on the differences between the CE and 
the GCE (in the context of ensemble averaged quantities) both due to quenched disorder and to thermodynamic 
fluctuations. The latter is manifest only at finite temperatures. This part includes careful examination of the notion 
of these ensembles in conjunction with the proper averaging procedures. Secondly, we shall address certain questions 
central to our understanding of statistical mechanics and thermodynamics on the mesoscopic scale. Thes questions 
include the following points: (1) We shall identify the (length scale dependent) temperature over which differences 
between the two ensembles are suppressed. (2) We shall relate this to the temperature (or length) scale over which 
differences between statistical mechanics and thermodynamics are suppressed. We shall thus identify the scales at 
which the thermodynamic limit is approached. (3) We shall relate the above to energy scales characterizing correlations 
of the underlying energy spectrum. These concepts will be illustrated below. 

The scope of the questions addressed here is quite general and includes generic many body systems. Our quantitative 
discussion, though, will be carried out for systems of noninteracting electrons. 

The outline of this paper is as follows. 

In Section 2 we discuss various definitions of "canonical" and "grandcanonical" . We distinguish between fluctuations 
in the number of particles (electrons), 6N, which arc due to thcrmodynamical effects, as opposed to quenched disorder 
induced SN. We also discuss the difference between canonical (or grandcanonical) constraints, defined with respect to 
the way the system is prepared (at equilibrium), and such constrains as defined when the system under consideration 
is subject to a dynamical perturbation. We next explain how to express ensemble averaged canonical quantities in 
terms of ensemble averaged grandcanonical quantities and fiuctuations thereof. This is done at zero temperature and 
is presented as an expansion to an arbitrary order in a small parameter. 

Section 3 contains a discussion of the differences between the CE and the GCE at finite temperature. This is related 
to differences between thermodynamics and statistical mechanics on the mesoscopic scale. Exact thermodynamic 
identities should be corrected when finite systems (at finite temperatures) are considered within the framework of 
statistical mechanics. This further modifies the relations between the CE and the GCE found above at T=0 alluded 
to above. We find how these differences between thermodynamic identities and statistical mechanical relations depend 
on temperature (and system's size) and evaluate the scale over which crossover to the "thermodynamic limit" takes 
place. It turns out that this scale (say, in temperature) is usually much larger than what one might normally guess 



(namely, that it is of the order of the level spacing) . Wc demonstrate our approach by evaluating the manner in which 
three different physical observables depend on the type of the ensemble employed; these observables are the heat 
capacity, the persistent current and the single level occupation probability. Differences between grandcanonically and 
canonically averaged quantities (e.g., the Fermi-Dirac distribution function and its canonical counterpart) are due to 
both quenched disorder and thermodynamic fluctuations. 

II. STATISTICAL ENSEMBLES AND QUENCHED DISORDER: T=0 

2a What Do We Mean by "Canonical" and "Grandcanonical" Conditions: Particle Fluctuations and 
Various Averaging Procedures 

The notion of ensemble averaging is in the very heart of statistical mechanics. Rather than considering a particular 
time evolution of a specific system with given initial conditions, we consider an ensemble of systems (or a set of initial 
conditions) and average over them with appropriate weights, stipulating that in any given set of experiments (or in 
an experiment that lasts a finite period of time) the experimental conditions are generally not reproducible to the last 
minute detail, but rather random sampling of the systems which form the ensemble takes place. 

The concept of a grandcanonical ensemble (as opposed to a canonical ensemble) invokes the notion of fluctuations 
in the number of particles, 6N. (In our case we shall refer to fluctuations in the numbe of electrons in the system). 
We now describe four different thought experiments which help us to eliminate some confusion that arises in the 
literature regarding what really is meant by CE or GCE. 
Experiment 1: We consider a system (in our context a mesoscopic conductor with disorder) coupled weakly to an 

external reservoir at temperature T and a chemical potential fi (Fig. la). The system will reach equilibrium with 
the external reservoir (if the coupling is weak the equilibration time may be long, but we are not interested here in 
transients). Even when equilibrium is attained the number of electrons in the system may still fluctuate in time; at 
finite temperature the system may exchange electrons with the reservoir. Denoting the number of electrons at times 
[3] ti,t2, ■ ■ ■ tn by A^i, iV2, . . . Nn we obtain a fluctuating sequence of numbers, whose variance is SN^^' . 
Experiment 2: We consider identical microscopic replica of the same mesoscopic systems, all weakly coupled to 

the same external reservoir (Fig. lb). Measuring the number of electrons in each system (at equilibrium) at a given 
arbitrary time we obtain a fluctuating sequence Ni,N2, ■ ■ ■ Nn whose variance is denoted as 5N^'^\ 

By the ergodic hypothesis SN^^'' = SN'-^K Moreover, they both tend to zero at T ^ 0. This reflects the fact that 
the origin of these fluctuations is dynamical: the Hamiltonian of the system contains a (small) coupling term to the 
reservoir. The number of electrons in the system is, therefore, not a conserved quantity. Evidently in the limit of 
T ^ 0, when the level broadening in the system due to the coupling term is smaller than the level spacing, and if 
fi is chosen to fall within the "gap" between two consecutive levels, the exchange of electrons with the reservoir is 
suppressed. 
Experiment 3: We now consider an ensemble of mesoscopic conductors which are macroscopically similar (they 

are made of the same material, are characterized by the same volume and impurity concentration, etc.) but are 
microscopically distinct (small variations in shape; different impurity configurations (Fig. Ic)). All these systems are 



weakly coupled, and equilibrated, with the same external reservoir. As a result, to each member of this ensemble 
corresponds a different spectrum (cf. Fig. 2a). Thus, even at T=0 (and with every member possessing the same value 
of chemical potential) different members of the ensemble will have different numbers of electrons. This defines SN^^' 
which docs not vanish even in the T ^ limit. The origin of this SN^^' is quenched disorder. 
Experiment 4: We may consider a single member of the ensemble at equilibrium, as function of an externally 

controlled parameter (magnetic field, Aharonov-Bohm (AB) fiux, electric field, external pressure, etc.) (Fig. Id). We 
shall denote this external parameter by <&. Generally the spectrum of the system depends on $. A generic situation is 
depicted in Fig. 2b. The number of particles (occupied levels) at equilibrium varies with $ (e.g., at T=0 A^(<i>i = 2, 
N{^2) = 3, cf. Fig. 2b), which defines SN'-'^K Again, this quantity does not vanish even in the limit T ^ 0. The 
order of magnitude of 5N^^' and SN^^' may be very different. For diffusive systems with $ representing an AB flux 
it has been found by Altshuler and Shklovskii [4,5] that for a two-dimensional systm (d=:2) and hard core scatterers 
SN'-^^ ~ h/EcT (for ( j)^ 3> 1) where t is the elastic mean free time and E^ = j^ is the Thouless energy {L is the 
system's linear size; i£ is the elastic mean file path, D is the diffusivity). By contrast [6] SN'-'^^ ^ 1. 

In the first part of this paper we shall set the temperature to be zero. We shall thus suppress any dynamic 
fiuctuations in SN, and the distinction between the GCE and the CE will be based on whether or not there are 
fiuctuations in the number of particles due to quenched disorder. (Only at a later stage shall we also consider finite 
temperatures). Thus, in the present context we shall employ the term CE to refer to a set of finite systems whose 
(sample specific) number of electrons N^ (i is used here as an index for the ensemble member) is selected independently 
of the disorder dependent energy spectrum, and is specifically $ independent. Similarly, a GCE implies that the 
sample specific chemical potential, /i*, is selected independently of the disorder energy spectrum, and is specifically 
$ independent. Throughout most of this paper we shall consider systems of noninteracting electrons, although some 
of the notions here may be generalized to include systems of interacting particles. 

It is possible to refer to a strong CE where N^ = N for all i, and similarly to a strong GCE where /i* = /i for all i 
(these were implicitly the conditions in Refs. 6-10 and Ref. 4, respectively). Alternatively, averaging within the CE 
(GCE) may be realized by randomly selecting values for N'^{^^) out of a broad and smooth distributions. We shall 
refer to these as weak CE (GCE) respectively. Below we shall argue that in many generic situations strong and weak 
averaging scenarios, as defined above, are equivalent, and will exploit this equivalence to our advantage. Figs. 3a, b 
depict schematically typical strong CE and strong GCE scenarios respectively. In Fig. 3a the probability that ji falls 
in a given interval between two levels is proportional to the distance ("gap") between these two levels. The effective 
(sample specific) chemical potential in the CE (T ^ A) is located half way between the last occupied (N—) level and 
the first vacant (A^ + 1—) level [11], Fig. 3b, regardless of the size of the interlevel gap. Thus, compared with the 
GCE, there is an enhanced probability to find narrow gaps in the vicinity of the effective chemical potential in the 
CE. The implications of this observation will become apparent below. 

The distinction among various types of ensembles may be carried further if we note that the nature of the ensemble 
depends conceptually on both the way it is prepared (brought into equilibrium), and the conditions under which the 
system is allowed to respond to an external bias. To be specific, let us consider an ensemble of multiply connected 
Aharonov-Bohm rings, subjected to magnetic fiux, <f>, threading each of them. At the preparation stage (which takes 



place at flux <I>), the rings may be brought into equilibrium either subject to canonical conditions {{N'^} are chosen to 
be random independent parameters, i.e., the number of electrons at each ring is selected independently of the ring's 
specific impurity configuration), or in accordance with grandcanonical conditions {fx^} are chosen independently of 
the impurity configuration). Subsequently, one performs measurements (e.g. of the ring's magnetic moments) at the 
fiux $. Then, too, as $ is varied (to obtain the response of the system), either canonical (number of electrons is 
kept unchanged) or grandcanonical (value of chemical potential is set fixed) conditions may be adopted. According 
to this classification, where a distinction is made between the preparation and the measurement stages, the CE is 
now denoted as canonical-canonical (C-C) and the GCE as GC-GC. We may also consider an intermediate, "hybrid" 
type of ensemble, GC-C. This "hybrid" procedure implies that the number of particles remains unchanged during the 
time $ is varied. This situation is achieved either when the system is decoupled from the reservoir before $ is varied, 
or if the coupling strength, 7, (affecting the inelastic broadening of the levels) [12] satisfies 7 < ujbias where w^^^ is 
the characteristic time scale for time modulations of the external bias. For 7 > ujbias the response to the external 
bias may be referred to as grandcanonical. We note that as long as we allow for (weak) coupling to the particle bath 
(which is necessary at the stage the system is prepared grandcanonically) and do not turn it off completely, the GC-C 
procedure is not an adiabetic one: when the rate of change of $ is slower than any other time scale in the system the 
procedure converges toward a GC-GC procedure. 

Following the above discussion and referring to the nature of both the preparation and measurement stages, we 
may consider GC-GC, C-C or GC-C procedures. At least in principle, one may realize these various ensembles 
experimentally [13,14]. The fourth combination, C-GC, is not physically relevant. Evidently the procedure we take at 
any given stage needs not be purely canonical or grandcanonical. For example, we may select the number of electrons 
within each ring to be only slightly correlated with the ring's disorder configuration. Thus there is a whole continuum 
of statistical procedures. 
2b Relations Among Various Statistical Ensembles: T=0 

Returning now to the two principle ensembles (CE and GCE, or, employing our alternative notation, C-C and 
GC-GC), it is desirable to find a relation between averages taken within the two ensembles respectively. From the 
calculational point of view, the GCE is more amenable to analytical studies (employing, e.g., disorder Green function 
techniques, random matrix theory (RMT), or supersymmetric (SUSY) calculations). 

To derive a general relation between GCE and CE averages we consider a strong GCE situation. The (thermody- 
namic average) number of electrons in the i^^ member of the ensemble is given by 

N\^,fi) = Jdei.\e,'^)f{e-fi), (1) 

where fi = const, and t^'(e, $) is a sample and flux dependent density of states (to be more concrete we presently 
assign to $ the meaning of an external magnetic flux). We define the ensemble averaged number of particles 

iV(M)^(A^^($,M)>- (2) 

The flux dependence of N{fi) is exponentially weak, hence we ignore it. We also define the mean level spacing as 



We further employ the approximation 



for A: > 2 (4) 



This is exact in two dimensions (d = 2) and involves corrections small in A//i in d ^ 2. 
For strong CE one has 

N^ jdev\e,'^)f(e-y}{<^)), (5) 

where N = const, and /i*(<&) is a sample and flux dependent chemical potential. Next we define the ensemble averaged 
chemical potential, p., as the solution of the following equation: 

N ^ N{J1) . (6) 

We then write the chemical potential as 

M*($)=^ + 5^^(<I>). (7) 

Expanding now the Fcrnii function in Eq. (5) and using the GCE relation, Eq. (1), one obtains 



oo ^ 



kl 

k=0 



d 
dp- 



l^N{^,p) 



iSf^-m)' . (8) 



Employing Eqs. (3), (4) and (6), we obtain the following transcendental equation for the dimensionless variable, 

OO 

X + ^ A'^X'^Cfc = , (9) 

where 



1 d'' 
kl 9/1* 



Ck = T^T;rj:SN; SN = N\^,p) - N{p) . (10) 



We shall look for the solution of Eq. (9) in the form of the series 

oo 

X = ^a„A". (11) 



n=0 



Substituting Eq. (10) into Eq. (9) and comparing respective orders of powers of A we find 



LJ_ (,„»_] , ,12, 



n\ 



where (n — 1) denotes the (n — 1)^* derivative with respect to /i. Finally one has a solution for (5/t'(<I>) 

(5MX'i>) ^ A y til:! A- (SN-^) ^""'^ . (13) 

This is a central result of our analysis so far, which enables us to calculate various canonical quantities. As an example 
we consider the situation where $ is an AB flux and work out the canonically averaged persistent current 
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n=0 

Employing a standard thermodynamic relation 

diGCE _ aiV'(^,M) 

and substituting Eq. (13) into Eq. (14) wc finally obtain 



^c.(<J>)-f:;Jy^^(V)". (14) 



d 






(15) 



(16) 



.ri=0 

One immediately notices that all terms in this series, except the n = 0, have the form 

a|0..)^ (17) 

This means that once averaged they yield a results proportional to A//!, as the only dependence of the averaged 
quantities for p, arises from the diffusion constant D — p,T/m (we are in d = 2, so the average density of states is 
constant). 

Equation (16) is a prototype of a relation connecting CE to a GCE average. It goes beyond the results of refs. 6-10 
in that it provides a systematic way to obtain higher corrections. More generally Eq. (16) can be cast in the form 

^\ ^^\ I A ^ V (-1)" .„ d^ 

d'^JcE 9^Jgce 9$^o(n + 2)! 9/2" 



^\ ^^\ +A—y ^~'> A-^SN-+^] (18) 



where $ now represents any static externally controlled parameter. Note that here SN is a grandcanonical quantity, 
so the canonical quantity on the l.h.s. of Eq. (18) is expressed solely by means of GCE quantities. 

A simple version of Eq. (18), namely this relation where only the n = term on the r.h.s. has been included, has 
been extensively employed [6-10] to study differences between the CE and the GCE persistent current in normal (as 
opposed to superconducting) rings and cylinders [15], and the anomalous paramagnetic orbital magnetism of small 
dots [16]. It is also interesting to consider the problem of "hybrid" ensembles, mentioned above [13,14,17]. 

III. STATISTICAL ENSEMBLES: T > 

So far we were interested in differences between the "canonical" and the "grandcanonical" ensembles (and hybrids 
thereof) due to quenched disorder statistical fluctuations. These differences are clearly manifested at T=0. The 
CE was represented as an ensemble of effectively grandcanonical systems, with sample specific (and flux dependent) 
chemical potentials. It is argued below that this is not the case at T > 0: even on the level of a specific system, and 
even when we consider thcrmodynamically averaged quantities (and not fluctuations), a canonical system cannot be 
replaced by an effective grandcanonical system. This will be the source to further differences between the CE and the 
GCE (i.e., differences beyond those discussed in Refs. 6-10,18). Here we study these difference and show that they 
are related to differences between statistical mechanics and thermodynamics on the mcsoscopic scale. 

In other words, consider an observable A. Canonical and grandcanonical averages will be denoted as {A^\n) and 
(A*|^) respectively. The difference between such averages can be written as 
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{A%) - (AX) = {A%^^^ - A^U) + {A% - A«|^.(A.)) . (19) 

The first term was the center of our discussion in Section 2. Note that if one assumes the relation 

^^U-^^U(A^,*) (20) 

(where $ is an externally controlled parameter) this implies that the second term on the r.h.s. of Eq. (19) vanishes. 
The first term on the r.h.s. exhausts, in that case, the difference between the CE and the GCE, including the 
temperature dependence. In general, though, as is discussed below, Eq. (20) is not satisfied exactly at T > 0, 
implying the non- vanishing of the second term on the r.h.s. of Eq. (19). The reason is that at finite temperature 
the relation given in Eq. (20) above, while being exact within the framework of thermodynamics, includes correction 
terms when considered within the framework of statistical mechanics and applied to finite systems (see Section 3b) . 
It follows that the fact that at finite temperatures we allow for dynamical particle fluctuations in the grandcanonical 
case (while no particle exchange with the reservoir is allowed under canonical constraints) affects quantities related 
to fiuctuation terms, e.g. {SN"^). In particular this applies to the relations between GCE and CE averages, which 
contain fiuctuation terms (cf. Eq. (18)). It is important to stress that the finite temperature difference between 
canonical and grandcanonical quantities, which is of a dynamical nature, exists even for sample specific quantities. 
This contribution is to be superimposed on top of the contribution discussed earlier, which arises due to averaging 
over quenched disorder. 

3a Differences Between the Canonical and the Grandcanonical Ensembles: Single Particle Level 
Occupancy 

We demonstrate the T > effect by studying the difference in the single particle level occupancy between canonical 
and grandcanonical systems. As a simple example, let us consider first a two-level system occupied by A^ = 1 electron 
(in the grandcanonical case this becomes N = 1). This is depicted in Fig. 4. Here A denotes the level spacing. At 
T==0 there is no practical difference between the canonical and the grandcanonical system. By contrast, at T > 0, 
the probability of finding the second level occupied is given by 

/c (level no. 2) — —^ canonical 

/gc (level no. 2) = . . grandcanonical (21) 

(in the canonical case we use the Boltzman factor for any given manybody state, projecting then onto the single 
particle state at hand; in the grandcanonical case we employ the Fcrmi-Dirac factor, / = foe- Clearly fc 7^ fee- 

We next consider a many level system. Closely related works [19] which include analytical studies of uniformly 
spaced spectra [20], possibly with the effect of level degeneracy [21], and a numerical study [22], have been published 
earlier. The Hamiltonian of the system, H, is assumed to preserve the number of particles (electrons), i.e. [H, N] ~ 0. 
The canonical partition function is given by 

Z^^^ (N) = Tr(e-^^'" S{N - N)) (22) 

where the hatted quantities are operators, and the Tr is unrestricted (i.e., it includes summation over eigenstates of 
N with various eigenvalues). We use the following representation of the Kroenecker delta function 



S(N-N)=p / i^e'3^(*-^) . (23) 

J 27ri 

Substituting Eq. (23) into Eq. (22) we obtain 

Zc{N)^(3 [ ^e-^^''TTe-^'("-'^^\ (24) 

Our goal now is to expand the r.h.s. of Eq. (24) around a saddle point [23-25]. To simplify the discussion we shall 
restrict ourselves here to systems of noninteracting fermions. We can find a diagonalizing basis in which H can be 
written as 



k 

The canonical partition function is then written as 



^ = E4''<«'^- (25) 



'Wi'm - R f _^ .-/3mA' TT /i ^e"'^('fc''~^) 






~3~ 

"'" givsiM) ^ (26) 



2m 



where 



g{y:) = -/3m + ^ E 1^"(1 + e-^('^''-^') • (27) 



TV 
fc 

In analogy with the grandcanonical ensemble, the single particle level occupation factor in the canonical case, 

fc\'^n;N) is given by 

f^c\el:^;N) = ^i-^Z«(iV) = -^-^^ . (28) 

At this stage we shall replace the limits of the integrals by infinite contours, /_'i^ -^ J_ , and evaluate these integrals 
employing the stationary phase approximation. It is possible to show [27] that this procedure is sound, and that the 
main contribution to the integral is contained within the original interval [^g^, ^] when the temperature is larger 
than the mean level spacing. 

One may now expand the numerator and the denominator of the r.h.s. of Eq. (28) (or Eq. (24)), making use of the 
diagrammatic expansion outlined further below. Here we shall derive a simple approximate expression, retaining only 
the Gaussian term in the exponent in Eq. (28), and expanding / up to the second order m^ — ^q (/Iq is the saddle 
point of the expansion, see Section 3b). One then obtains (on the limit of a quasi continuous spectrum): 

/c(e; TV) = /{e m^^^) - ^ /"(^ - M^^^) • (29) 



We readily see that the canonical distribution fc is qualitatively equivalent to the grandcanonical distribution at a 
lower temperature. To find this effective lower temperature we write 



This yields 



f'i^l■, Pgc) = /c(e = m; N; (3c) = /'(m; Pc) ~ ^/"(a^; Pc) ■ (30) 



/3gc«/3c(1 + ^) , (31) 



where (3qc is the effective grandcanonical inverse temperature which corresponds to a canonical inverse temperature 
(3. Based on the above discussion, the last expression is valid for T ~ A. It can be rewritten as 

Tgc « ^^V (32) 



which in the limit of Tc 3> A becomes 



1 I 4Tc 



Tgc~Tc-- Tc:^A. (33) 



One may employ slightly different criteria for defining the effective temperature of the canonical distribution function 
(e.g., /(/i,/3Gc) = fci^ — M;-^;/^c) instead of Eq. (30)). But the statement that the canonical distribution, at 
temperatures larger than A, corresponds to Tc > Tec should remain valid. A similar shift of the effective temperature 
has been discussed (in the context of a uniform spectrum) in Rcfs. [20-22]. This also agrees qualitatively with the 
2-level scenario discussed above, where the effective canonical temperature is higher by a factor of 2. 

3b Differences Between the Canonical and the Grandcanonical Ensembles at T > 0: A General 
Approach 

Our next step is to develop a systematic method for calculating physical observables within the canonical ensemble, 
accounting for finite temperature effects, in particular for the fact that at finite temperature fc =/= ,f ■ We incorporate 
into Eq. (24) for the canonical partition functin, Zc, the relation 

writing 

Zc(iV,T,<&)«^ / d^e-^(^(^'^'*)+^^) , (35) 

2m J 

The limits of the integral have been extended to infinity, following the same arguments as in the derivation of fc- 
We note that the expression for Zc is quite general and is valid for systems of interacting electrons as well. We shall 
restrict ourselves to systems of non-interacting electrons in the diffusive regime. Our goal here is to evaluate ensemble 
averaged thermodynamic observables within the canonical ensemble. (We strec that similar ideas may be employed 
for calculating other classes of observables, e.g., linear response transport coefficients). 

As we have done above, we define the sample specific saddle point of the integral in Eq. (35), ji^ [T, $) [28] 
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dy. 



[0«(Ai,T,$)+^7V] 



= 0. 



A'^^'CT,*) 



Let us define the grandcanonical expectation value of the particle number as 



M^'(.,T,.)^-^~il. 



(36) 



(37) 



(For discussion of the relation between this identity, defined within statistical mechanics, and thermodynamic identi- 
ties, we refer the reader to Section 3d). The saddle point equation then assumes the form: 



iv«(M^''(T,$),r,$)=iv, 



(38) 



where A'^ is a given integer. We are supposed to solve Eq. (38) for /ig , and substitute it in the expression for Zq, 
Eq. (35), expanding around the saddle point: 



r!«(M)+M^-F«(Mo'') + 



1 92r!( 



2! 9^2 



{^.-^^'Y + ...+ 



1 9"f7(*) 



{p-f^^o^r + .. 



(39) 



We note that /ig is real. We have used the notation F^^^{iiq ) = f7'^*'(/io)+/io jV. Unfortunately the above expression, 

Eq. (39) is not very useful for practical purposes. Each term is to be evaluated at the chemical potential /Iq which 

is a function of the external parameter $, but even more significantly, is sample specific. Evidently, it is important 

to perform ensemble averaging by calculating the sample specific function /ij, ($) for each and every member of the 

ensemble. Instead we replace this quantity by a (positive) constant jl, assumed to be the result of averaging fi^ ($) 

— : * 

over $ and over ensemble realizations, fl = (/ig ($) ). We assume that in some sense (to be discussed below) p, is 

close to /Xq ($). We thus distort the integration contour, C, of Eq. (35) which goes through the saddle point /ig ($) 

into another contour, C, which goes through fl. Expanding around the latter, Eq. (39) is now replaced by 



f^« (m) +fiN = F« if,)--Li^-fif-J2 



dp" 



-SN'^'^ 



(m-m)" 



(40) 



Here SN'-^^ = A^(*)(/2) — N. The mean level spacing A satisfies ^ — { ^1 - ) -In writing Eq. (40) we assumed that 
N — -^j such that -§^^N = for n > 2. Strictly speaking this is correct for a two dimensional system; for d 7^ 2 we 
introduce errors of order ■= or less. 

Substituting the expansion, Eq. (40), into the expression for Zc, Eq. (35), we obtain 



Z«(7V,T,$) = ^e-''^'''(^)|d,.e^(^-^) .exp|/3f;l^^^^(M-M)" 



2m 



(41) 



We now define new variables 



-*a/^(a*-m) 






(zY 



n\ (/3A)t-i 5^ 

in terms of which we may rewrite the partition function as 



nn-l 

A"-i4^--j-,57V«(/i) , 



(42) 
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4H^>r,$).y^e-^-'*'(^) /-^e^ expjfl/ivj 



(43) 



Equation (43) is the basis of our diagrammatic expansion. We first stress that ah diagrams in this expansion are to be 
evaluated within the grandcanonical ensemble, which is the attractive feature of our analysis. The vertex Vn should 
be understood in the following way: SN is a trace of a ring diagram with one scalar vertex (i.e. proportional to the 
first derivative of fl with respect to p.). Hence Vn corresponds to n scalar vertices. (One may dress such vertices by 
interaction lines, impurity lines, etc. as depicted in Fig. 5). We now expand the exponential in Eq. (43). We shall 
comment on the small parameter of the expansion below. In analogy with Wick's theorem we are now able to define 
contraction with the variable t. We shall employ the following notation: for the second order term in r we write 

This defines the combinational factor that comes with this term (=1), corresponding to the fact that there is only 
one way of contraction. For the fourth order term in t we write 

dT 



'2. -^^'- <«' 

The combinational factor associated with this term is 3, corresponding to 3 possibilities of contraction. In these 
"contractions" the role of the zeroth order propagation is played by e 2 . In our diagrammatic representation we 
shall use a zig-zag line for this propagator. We stress that the only role of this statistical r-propator is to account for 
combinatorial factors. It does not carry any momentum or energy. We also note that odd powers of r vanish. 

To evaluate thermodynamic derivatives we need to calculate log Zq. To this end we may employ the linked cluster 
expression (including statistical lines!) whereby only linked diagrams are to be accounted for (cf. e.g. Ref. 29). We 
thus may write 

i^(*)(Ar,r,$) = F!^%(fi,T,<i>) - -^log\ — - {all connected diagramsY"^ . (46) 

p V 27r 

Anticipating the discussion of the next section, partaining to differences between statistical mechanical and thermo- 
dynamic quantities, we use the notation 

4^(A2,T,$) = r!«(^,T,$)+M7V. (47) 

The first term on the r.h.s. of Eq. (46) yields the sample specific grandcanonical contribution. The second term, due 
to Gaussian fluctuations around the expansion point (not the saddle point!) /I, is a temperature-dependent sample- 
independent contribution. All the other terms are, in principle, sample specific. The quantity r2(*-'(/i,T, $) may be 
evaluated employing the standard linked cluster expansion [30], which does not include statistical lines. The r.h.s. of 
Eq. (46) yields terms which are particular to the canonical ensemble. 

At this point we can formulate the rules for constructing canonical diagrams. We consider a linked diagram 
consisting of p electronic loops (with ni, n2, . . . rip scalar vertices respectively). The total number of vertices must be 
even ni+n2 + ■ ■ -Up — 2K. As all these vertices are connected by statistical lines, it follows that K is the number of 
such lines, K > p — 1. 
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1. With each statistical hne associate a factor (—A). 

2. With each electronic loop which contains Ui vertices associate a factor 

3. Any subset of m loops which consist of the same number of 
vertices carries a factor ^. 

ml 

4. Each diagram carries a factor T™, w ^ K — p + 1. 

5. With each diagram associate a combinatorial factor, representing the 

number of different ways of constructing it. (48) 

We stress that the numerical factors obtained through the contractions discussed above correspond, in the present 
language, to the number of different ways of connecting scalar vertices by statistical lines. Hence, the only role of the 
statistical lines is to provide the right counting of diagrams. We also stress that the tempertaure factor T"" associated 
with each diagram (see rule no. 4 above) does not exhaust the temperature dependence of that particular diagram. 
The latter usually contains further important temperature dependence. We shall nevertheless use this power of T 
to classify the diagrams into different families, characterized by the index w. Examples are depicted in Figs. 6 and 
7. It is understood that disorder averaging has to be carried out subsequently. At zero temperature only the w = 
family survives. It is easy to see that these diagrams originate from the linear term in Eq. (40). An expansion about 
the true saddle point, /ij, (T, $) would not reproduce these diagrams. It turns out that these are the contributions 
which describe differences between the CE and the GCE at T = 0, cf. Section 2. They are obtained when Eq. (20) 
is assumed to hold, and each member of the canonical ensemble is assigned an effective, sample specific chemical 
potential. The w — family represents contributions due to quenched disorder (as in Refs. [6-10,14]) but not due 
to differences between statistical mechanics and thermodynamics (see Section 3d). After some algebra one obtains 
(hereafter we suppress the index (i) and consider only ensemble averaged quantities) 

n—l ^ ' 

Derivation of Eq. (59) with respect to an external parameter, <I>, yields Eq. (18) [14]. The first (n = 1) term 
in the sum corresponds to the two loop diagram in Fig. 6. Upon averaging over (diffusive) disorder, it yields the 
Altshuler-Shklovskii term [4], which was employed in Ref. [6-10]. The n > 2 terms include complete derivatives with 
respect to p,, and are negligible upon averaging (being small in the parameter A//i). This provides us with a posteriori 
justification of the diagrammatic expansion: sample specific terms in this expansion are not necessarily small, but the 
ensemble average is well-behaved. 

We next include the w > 1 families too. It can be shown that for the regime where we employ our expansion, 
T > A, the leading term (in A/T) of each family is represented by a two-loop diagram (cf. Figs. 6,7). Evaluation of 
these diagrams (neglecting full derivatives with respect to p) results in 

A^ (AT)- /fd-'SNW 
(<5F,.<,_,<,„,) ^-^ ^-^^ (^ j j . (50) 
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The w = term is the Altshuler-Shklovskii quenched disorder contribution [4]. 

3c Relation to Spectral Correlations: Persistent Current and Heat Capacity 

The summation in Eq. (50) is clearly related to correlations in the spectra of finite size disordered systems. We 
define the correlator /C and its Fourier transform JC: 

/C(e - e') ^ A2((z.(e)Ke'))) - (He)) {{e'))) ^^ J rfi^(t)e^*(-^')/^ , (51) 

— oo 

where i>{e) is the sample specific density of states and t is the dimensionless time (in units of h/A). We can write 



„ +00 

8^i)N\ \ 8^ 8^ f [' dfdf' 



(52) 



where / is the Fermi-Dirac function. Changing variables to ^ = e — e', ry = (e + e')/2j performing the integral over 77, 
and Fourier transforming with respect to ^, we obtain 

W_.„,) .-^ [-) j^t-^^^{ij) . (53) 



W = g 



Equations (46) and (53) form the basis for the analysis of the various corrections to the GCE averages, and depict 
the dependence of these corrections on spectral correlations. We study two examples comparing the w — (quenched 
disorder contribution) and w — I. We consider the case of diffusive disorder. 

The average persistent current in the canonical ensemble, {I)ce, is obtained by deriving (F) with respect to the 
Aharonov-Bohm flux, $. The flux dependent part of the time correlator for quasi-one-dimensiional rings is given by 



1,1 / -1 U<J 2 

where g — Ec/A. 3> 1 is the dimenionless conductance (Ec is the Thouless correlation energy) and $0 — he/ e. Only 
even harmonics appear in Eq. (54). Expanding {I)ce = SS^i-^pSin47rp(f>/$o, we find that the quenched disorder 
contribution (w — 0) is [6-10] 

We shall denote the leading w — 1 contribution with a superscript SM — TD. The rationale behind that will become 
evident in the next section. We obtain 

The average canonical heat eapaeity {{C)ce)- It is given by {C)ce — ^0^d'^{P{F))/d0^ and is written as a sum of 
contributions (cf. Eq. (46): 

{C)cE = {C)gce + {SC'^^^'n + {SCn + {6C^^-^'') . (57) 
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The first term is tlie grandcanonical contribution, which for a degenerate gas of non- interacting electrons is ^ -^ . The 
term due to Gaussian fluctuations yields —2- This contribution can be reinterpreted when T > A as a shift of the 
grandcanonical temperature (T > A) towards a lower temperature T ^ T — 2^2- A forced by the canonical constraints. 
This is qualitatively similar to the shift found in Section 3a, in the context of the single level occupation function. 
A similar shift in the heat capacity was found in Ref. [20] (see also Refs. [21,22]). The next terms in Eq. (57) are 
evaluated employing Eq. (53). We are interested in energies larger than A, hence t < 1. The time correlator has two 
interesting regimes (random matrix theory and Altshuler-Shklovskii [4] respectively): 

where b = 1,2,4 for the orthogonal, unitary and symplectic ensembles respectively. Here t^i is the elastic mean free 
time. Eq. (58) does not account for the crossover regimes. It leads to 

and in the leading order in A/T [31] 

2 ( KM A <T < E 

In this example the SM — TD contribution is parametrically smaller than the contribution due to quenched disorder. 

3d Differences Between Statistical Mechanis and Thermodynamics 

Previous studies of differences between the GCE and the CE used the thermodynamic identity, Eq. (20) as a starting 
point. This implies that for any specific realization, a canonical system can be replaced by an equivalent grandcanonical 
system with an effective chemical potential which depends on disorder and on the value of the relevant external 
parameters ($). Physical observables, though, should be evaluated within the framework of statistical mechanics. 
While Eq. (20), considered as a thermodynamic relation, is exact - by construction, a statistical mechanics relation 
between F and fi contains important corrections, cf. Eq. (46). These corrections are expected to be suppressed in 
the "thermodynamic limit" - at high enough temperatures and/or large enough system's size. 

These types of corrections (deviations from the thermodynamic relation, Eq. (20)) may be quite important on the 
mesoscopic scale. Thus, the expansion outlined above includes two qualitatively different classes of contributions: (a) 
the w — Q family, which assumes the validity of Eq. (20), and accounts for differences between the CE and the GCE 
due to quenched disorder; (6) the w >1 families, which account for deviations (within statistical mechanics) from the 
thermodynamic identities. (This is the rationale behind the notation used in Eqs. (56), (57) and (60)). The two-loop 
term of the w = 1 family is the leading contribution (in A/T) of this class. 

Coming now to the examples discussed in the previous section, we note that the ^^SM—TD^^ contributions arc smaller 
than those due to quenched disorder. It is interesting to note, though, that the former is quite robust as function 
termpcrature: contrary to what might be a naive expectation, contributions due to differences between statistical 
mechanics are not suppressed on the scale T ^ A. In the case of the persistent current the "5M — TD^^ contribution 
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(Eq. (56)) decreases algebraically with temperature up to T ^ Ec ^ A, and is then suppressed exponentially. In the 
case of the heat capacity it is also suppressed algebraically (up to T ^ fi/Tei- Our analysis underlines the relation 
between the ''SM — TD" contributions and spectral correlations. 
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Figure Captions 
Fig. 1: Four different thought experiments leading to different definitions of 5N (see text). 

Fig. 2: Fluctuations in the number of electrons done to quenched disorder, (a) SN^^': At zero temperature members 
(1), (2) and (3) of the ensemble (differing in their respective shapes, impurity configurations, etc.) will have 3, 
4 and 2 electrons, respectively, (b) SN'-^^: Fluctuations in A^ due to variations of an external parameter, $. 
The equilibrium values of iV at $1, $2 are 3,2, respectively. 

Fig. 3: Empty-occupied levels near the Fermi energy, (a) GCE; (b) CE. Enhanced probability for chemical potential 
to fall in a narrow gap. 

Fig. 4: A two-level system under (a) canonical, (b) grandcanonical constraints. There are no observable differences 
at T = 0. 

Fig. 5: A contribution to V5. Solid lines denote electron propagators; wiggly lines - electron-electron interactions; 
dashed lines - impurity scattering. 

Fig. 6: The w — family of skeleton canonical diagrams with K = 1,2,3. Zigzag lines are "statistical lines"; full 
lines - electron propagators; black dots - scalar vertices. The first diagram yields the Altshulcr-Shklovskii result 
upon averaging over disorder. 

Fig. 7: Skeleton diagrams of the w = 1 family (up to K = 3). 
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